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Abstract. Two recent investigations are reviewed: quantum effects for DNA aggregates and 
scars formation on virus capsids. The possibility that scars could explain certain data recently 
obtained by Sundquist's group in electron cryotomography of immature HIV-1 virions is also 
briefly addressed. Furthermore, a bottom-up reflection is presented on the need to invent new 
physics to pave the way to a rigorous physical theory of biological phenomena. Our experience 
in the two researches presented here and our personal interpretation of Schrodinger's vision are 
behind the latter request. 



1. Introduction 

It seems to us that the frontier of physics is rapidly moving to the biological side not 
only in the obvious way that the highly detailed knowledge of the biological data requires 
increasingly sophisticated physics for their explanation, but also (and even more intriguingly 
for the theoretical physicist's mind) in the sense that the time might be mature for a profound 
reformulation of biology that makes it a science similar to physics: few principles, mathematical 
derivations, explanation and prediction of phenomena (and back). Nonetheless, the problem is 
very difficult if faced with the proper rigor and we have no answer to propose, a part from the 
conviction that some important changes are necessary on the physical approach to the problems 
of biology and that a "physics tailored for life phenomena" needs to be invented. On this we 
comment only in the last section where we arrive by taking a feet-on-the-ground road that makes 
two stops prior to getting there: the study of (relativistic) quantum effects in DNA condensation 
P] and the role of certain topological defects (scars, see, e.g., [2]) for virus capsids shape-changes 
[3]. This way we have the chance to face important biological problems on the table and to try 
to learn on the battle-field what is that is missing from the big picture. 



2. Quantum fields holding together DNA aggregates 

2.1. The story 

DNA molecules in aqueous solution ionize and become highly charged anions which strongly 
repel |3|-[l2] (for reviews see, e.g., [THE]). When specific cations are added (i) DNA attracts 
and binds them to make a new structure, the DNA-cations complex, and (ii) when about 90 
per cent of the DNA negative charge is screened the like-sign DNA-cations strands attract and 
collapse (see, e.g., [E2 EH ES]) to form finite-size aggregates whose shape is either rod-like or 
spheroidal or, most commonly, toroidal |16| . 

Long ago Oosawa and Manning (OM) explained the counterions condensation as a phase 
transition within the classical Poisson-Boltzmann (PB) mean field electrostatic theory |17j . 



Despite important advancements, the second part of the puzzle still presents many open 
questions: (i) it is not understood why the aggregate does not grow forever; (ii) there is no 
general consensus on the necessity to go beyond the PB classical theory, as proposed in [18] 
(see also [5]) and as opposed in [7J (see also [HI [19]); (iii) whichever approach is used (the zero 
frequency Casimir/van der Waals interaction [20J (see also }21|); the "electrostatic zipper" model 
of [TJ; correlations of thermal Gaussian fluctuations of the number density of counterions [5]I18|: 
the Wigner crystal approach [llj [22]) the paradigm is that the interactions are classical |13| . 
The lack of appreciation of quantum effects for this phenomenon was probably due to the lack of 
calculations based on the appropriate codimension two (lines in three dimensions) Casimir-like 
techniqu^l] [25j and with cations not directly participating in the interaction but playing the 
passive role of screening the electrostatic repulsion. 



2.2. A new model based on relativistic quantum fields 

In pp we singled-out the zero-point interaction due to the disturbances induced in the quantum 
electric vacuum by the presence of the (nearly) neutral DNA-cations complexes. We modeled the 
N anions (the DNA strands) as infinite lines (finite length effects are unimportant) all parallel to 
the z-axis and located at Zj in the x — y plane with certain coefficients Viizi) carrying information 
on the charge structure of the DNA and taken to be Vi(zi) = v = constant, Vi = 1, ...,N. The 



cations screen the total charge and set the length scale, 
potential <£(x) in this set-up is given by [1] 
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which is a modified Debye-Hiickel (DH) equation, with fj? = k 2 K 2 , k the cations valency, 
= (eksT/ (87re 2 no)) 1 ^ 2 the Debye screening length, e the dielectric constant of water, no fixes 
the zero of the potential, A = 4irv\q\ 2 /ek B T, J = -(l/e)4n\q\u J2?=i ^ (2) (%± ~ k), T = 300K, 
and with Eq.(pQ) the limit for weak $ of a PB equation modified to include the DNA strands in 
the charge (Boltzmann) distribution [1] 
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where q < is the charge of the DNA strand with n,£) NA (x) = J2iLi Vi(zi)5 2 (x± — k) and N the 
number of strands. 

Eq.(pQ) is only a starting point, as our concern are time dependent fluctuations, ^(x) — > 
&(x) + 4>(x, t), that we treat as quantum by considering the associated action (we set % = c = 1, 
c speed of light in the medium) 
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which we infer from the action for $ that gives Eq.([T|) as equation of motion, A(<&) 
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A Eili o~ {2) (x± - li)}$> + J$) , to be of the order of h. Note that in 



A(<p) the term with the coupling to J is zero because / d 4 xJ(p = J d?xJ Jq dt(fi = 0. 

These fluctuations are disturbances of the quantum vacuum induced by the delta- 
functions and, since they are the electric field, they travel at the speed of light in the 



1 The work in |20l [6] is based on the Lifshitz computation for the Casimir effect in presence of dielectric media 
|23] (see also [24] and references therein) i.e. a codimension one (surfaces in three dimensions) calculation. 




Figure 1. Two-Body interaction energy /length £{x) against lattice size x for a = 1 (upper 
curve) and a = 2 (lower curve) 



medium. We considered the generating functional Z[<&, <^>] = flD&je 1 ^^ J[D(f)]e ij ^^' = 
/[D < £]e~(' /l (*) +correctlons ) where we Wick rotate on the time direction t — ► it, and by identifying 
the effective action as «4eff(5>) = A($) + £t one sees that pQ 

where /o(-E') is the density of states, i.e. £ is the zero point "Casimir-like" energy of the system, 
that is found to be p], [25] (see also [26J) 
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where we reintroduced fi and c, Kq(x) is the modified Bessel function of the second kind of order 
zero, \x is the scale parameter introduced earlier, Zj,- = |Zj — L| are the relative distances and 
M is a further mass scale parameter originated by the codimension two (for stability it satisfies 
M < //). We are not considering charge density or positional fluctuations of counterions, hence 
we do not have slow-moving fluctuations, as, e.g., the authors of [6] that had to consider only zero 
frequency (i.e. classical) contribution to the usual Casimir/van der Waals effect. We have here a 
universal mechanism of interaction that (i) only becomes important when the OM condensation 
has taken place, (ii) whose temperature dependence is via the length scale // and (iii) with a free 
(but highly constrained in units of fi) parameter M. 
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2.3. New results 

The two-body interaction indeed is attractive as shown in Fig. [TJ Our focus needs be on 
configurations as close as possible to real cases [16], such as that shown in Fig. [21 where the 
relative distances are Z« = CijX with x the basic lattice distance and Cj,- taking the symmetry 
of the given arrangement into account. The units are fi^ 1 ~ O(10) A, for distances, and 
£ 10" 2 eV/A, for lineal energy density. M = e l l a (we choose a = 1 and a = 2). Thermal 

fluctuations, as computed, e.g., in [5], give for the two-body interaction a maximum value of 
£ ~ 5 x 1CP 3 eV/A at 10 A. At this distance our two-body interaction is 5 x 10" 2 eV/A (for 
a = 1) or 2 x 10 _1 eV/A (for a = 2), i.e. between one and two orders of magnitude stronger. For 
the many-body case, the case of importance for the aggregates, this factor grows enormously [1] 
but we cannot trust our approximations for x too close to the singular value of the logarithm. 
It is clear, though, that at the distances of relevance this quantum energy is stronger (or much 
stronger) than thermal energy. Thus we can conclude that it is relativistic quantum field zero- 
point energy that holds together DNA aggregates! This is nice but... 




Figure 2. Configuration used to compute the interaction of 19 strands. 
2-4- Effects of the lack of a general theory of biological physics 

... important pieces are missing here due to the lack of a proper "general frame": (i) where is, 
in the physics of the problem, that those filaments are the carriers of genetic information? (ii) is 
the description of this phenomenon complicated because we are not using the most appropriate 
approach? and is that meaning that we need new variables or a different coarse-graining or that 
we need to take into account that DNA compaction is just part of a series of "moves" made with 
a "purpose" (e.g., to pack the DNA in the bacteriophage's head or into the sperm cell and that 
is, in turn, just one step of virus assembly or of fecundation, and the latter one step in the "life" 
cycle of the virus or in the reproduction of eukaryotic cell organisms, and so on to higher and 
higher levels of complexity) (iii) how can we translate into physical language the "purpose", or, 
in other words, the clear fact that all these moves are done by "using the laws of physics for a 
goal" ? is it matter of changing the definition of physical rules or what? finally, (iv) is quantum 
physics used to full strength, i.e. as pivotal for life itself [27] (see also [2"8])? 

3. Scars on viruses: the general conjecture and the case of immature HIV-1 

3.1. The two stories: virus structure and Thomson problem 

A virus is a piece of DNA or RNA surrounded by a protein coat, the capsid, sometimes cased into 
a lipidic membrane, the envelop [29J. Capsid's shapes can be: helical, icosahedral or complex 
(sphero-cylindrical, conical, tubular or more complicated shapes) [29]. Viruses may change their 
shape (polymorphism) as an important step in their "life"-cycle. For instance, HIV-1 is only 
infective when its capsid has changed from spherical to conical (maturation) [30J. 

The theory of icosahedral capsids was proposed by Crick and Watson [31] and later established 
by Caspar and Klug (CK) [32]: those capsids all have 60T proteins arranged into 12 pentamers 
and 10(T — 1) hexamers, where T = n 2 + m? + nm = 1, 3, 4, 7, .... These numbers descend from 
the Poincare-Euler theorem!! 

J2(6-P)N P = 6 X , (6) 

N 

applied to the sphere, x = 2, and from the CK "quasi-equivalence" principle [32], see 
also [3J. The CK theory is nowadays an established paradigm among virologists [5U] and 
various modifications/generalizations have been proposed by physicists and mathematicians 
[21 [331 EH EHl [36] (for a review see [37]), but none made its way to virology textbooks. 

An intimately related physical set-up is that of N electrons on the surface of a sphere. To 
find their minimum energy configurations means to solve the so-called "Thomson problem" 

2 N p is the number of p-gons used to tile a surface, e.g., N$ pentagons, N§ hexagons, etc. The resulting polyhedron 
P has Vp — 1/3^2 N N p p vertices, Ep = 1/2^2 N N p p edges, and Fp = ^2 N N P faces, giving for the Euler 
characteristic \ = Vp — Ep + Fp the expression in Eq. (|6]). 
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Figure 3. Idealized scar formation-annihilation mechanism. 




|38j . unsolved in general. For N < 0(500) (and of the form N = 10T + 2) the icosahedron 

is the solution, as shown in [39], a work inspired by the CK theory of virus capsids. For 
bigger N configurations that present "scars" are the solution as shown in [40J and for spherical 
elastic materials (non-Coulomb potentials) in [IT]. A scar is a lineal pattern of pentagons- 
heptagons that initiate at the pentagonal vertices of the icosahedron, whose length (i.e. number 
of pentagons and heptagons in it) and geometrical arrangement depend on N (see [2j for an 
approach based on spontaneous symmetry breaking of the icosahedral group). Scars have 
been experimentally found (in spherical crystals of mutually repelling polystyrene beads self- 
assembled on water droplets in oil [12]) but still need a deep understanding in terms of phase 
transitions. Nonetheless, it is easy to see that they are allowed by the theorem ([6]) written as 



This means that Nq can be arbitrary (hence also the required N$ = 10(T — 1)) and N^ — Nj = 12, 
hence, starting off with the 12 pentagons of the icosahedron one can add pairs pentagon- 
heptagon, but not a pentagon or a heptagon separately. Geometrically, we are saying that 
a unit sphere has curvature i? sp here = +1 and each polygon contributes to this curvature with 
R p = (6 — p)/l2: a hexagon with Rq = 0, a pentagon with R§ = +1/12, a heptagon with 
i?7 = —1/12. The total energy is then the sum of bending energy and stretching energy, 
£ t = £ b + £ s that compete for the minimization of the total energy [41j . 

3.2. Scar formation- annihilation mechanism for viruses 

Below the threshold for the scar production N scaT , £b is provided by the 12 pentagons, while the 
hexagons contribute to £ s only. One way to imagine the transition at N > N BcaT is to think of 
a pair 6-6, with zero total and local curvature and zero bending energy, converted into a pair 
5-7 (see Fig. [3|), again with zero total curvature but with nonzero local curvature hence with 
nonzero bending energy given by 2Eb, where is the energy necessary to convert a 6 into a 5 
or into a 7. In [3] we proposed that scars can appear on virus capsids at an intermediate stage of 
their evolution towards a non-spherical shape and can actually drive such shape-change. There 
Eb should be related to the conformational switching energy [33] (for HIV-1 see, e.g., [S]). 

The (idealized) shape-change mechanism of [3] is (see Fig. [3]): i) proteins first make a CK 
icosahedron; ii) when, e.g., N reaches N SCSuT , they form scars: 6 — 6 — > 5 — 7; iii) the capsid 
changes shape via the release of the bending energy into stretching energy at the location of the 
scar with the consequent "annihilation" of the 5-7 pair: 5 — 7 — > 6 — 6. The resulting capsid has 
the CK morphological units but not the spherical shape. 

An example given in [3] is that of scars created only near the 10 inner vertices via a mechanism 
that respects a C5 rotation symmetry around the north pole-south pole axis. In Fig.H]the vertices 



(6 - 5) N 5 + (6 - 6) N 6 + (6 - 7) N 7 = 12 . 



(7) 



where the scars are formed are indicated with •, while the other two are indicated with o. If 
we require that this mechanism is area preserving, i.e. the total number of proteins before and 
after is always 60T the final capsid obtained is the spherocylinder of Fig. 5 (a shape taken by 
certain bacteriophages' head) with the 12 pentamers at the vertices and the 10(T — 1) hexamers 
distributed differently with respect to the intermediate icosadeltahedron. In general a variety of 
final capsid shapes could be obtained via this mechanism. For the C5 symmetric example given, 
if, for instance, the scars carry a bigger bending energy the final shape could be a backy-tube 
capsid; if the orientation of the scars in the previous setting is such that C shrinks, hence B 
becomes longer (see Fig. 5) then a disk-like shape is obtained; etc.. 




Figure 4. The intermediate Figure 5. The final spherocylin- 

spherical capsid. drical capsid. 



3.3. Are scars at work on the HIV-1 immature virions? The data of Sundquist's group. 
This is just an idealized model and to make full contact with real viruses one has to be ready 
to adjust it and, as scars are supposed to disappear, it is also a delicate matter to find the 
right experimental set-ups for spotting the scars. A case that seems particularly promising is 
the important case (for human health) of the HIV-1. There the number of Gag polyproteins 
(essentially made of the three proteins MA, CA, NC) that make up the immature quasi-spherical 
capsid is estimated to be between 2500 [35] and 5000 [36] (even 11000 is reported [46]), i.e. well 
above the iV scar of the Thomson problem. Indeed, the impressive data of Sundquist's group taken 
by electron cryotomography of the immature capsid [35] (see also [30]) reveal the structure of 
the CA shell that is seen to have an hexameric texture and large areas with an apparent lack of 
symmetry ( "disordered " areas in the terminology of [35] ) . The role of these a-symmetric areas in 
the maturation process, the apparent lack of pentameric structures, the mechanism that allows 
for closure of the spherical capsid or the absence of such closure, and other aspects are important 
issues to be understood. 

We claim here that in those areas of the capsid of the immature HIV-1 virions the scar 
formation-annihilation mechanism we discussed might be at work, although, since it takes place 
within a many-layer structure comprising the rest of the Gag and the lipidic envelop, it does not 
give the effect described above in the idealized case. What it presumably does instead is: (i) to 
stretch the structure, hence not allowing the perfect spherical shape; (ii) in combination with 
the proteolytic cleavage is responsible for the dramatic reduction of the number of proteins (the 
generally accepted number for the mature capsid is around 1500); (iii) it might even explain the 
apparent absence of pentamers if one supposes that when the annihilation takes places it involves 



also the "primitive" pentamers, i.e. what in the model are the untouched (by the annihilation) 
vertices of the intermediate icosahedron. 

Of course we need to explore those assertions in a quantitative fashion within the model, and 
we have work in progress in that direction. In the meanwhile, a faster way is to probe those 
claims directly in experiments, ideally with the help of Sundquist's laboratory. 

S.J,.. Effects of the lack of a general theory of biological physics 

Again in this case, as for the DNA compaction, it is clear that the appropriate general frame is 
lacking. With some adjustment (e.g., "DNA packing" goes to "capsid shape-changes", etc) all 
questions of Subsection 12.41 apply here, except (i). One important point here is that viruses are 
at the border between dead matter (where standard physics applies) and living matter (where 
it appears that new physics is necessary), hence we deem them to be the best candidates to 
understand how to build-up the missing frame. For instance, invariably all models (including 
ours) on virus structure deal with (free) energy minimization (see [17] for a review), just like, 
for instance, an elastic ball. May be this is a place to start looking for changes. 

4. Biology as Physics 

So, what is it this new physics that we claim to be necessary for a true understanding of living 
matter? The vast majority among physicists and biologists says that there is nothing to be done: 
the physics for living matter is just the same as for dead matter (elastic balls, charged wires, 
classical machines at the nanometer scale, etc.) it is "just more complicated to apply". Even 
taking this conservative view, it seems to us that (to say it with a metaphor) it is like trying 
to solve a (difficult) geometry problem with spherical symmetry insisting in using cartesian 
coordinates: results are more difficult to get, their interpretation becomes obscure and, most 
importantly, we might miss global information insisting with a description that is good locally 
but not globally (see also the metaphor of the mechanical engineer bumping into an electric motor 
of [21] )■ The conservative approach to biophysics is nowadays achieving impressive results, and 
this approach is just what we have been using in our own work reviewed in the previous two 
Sections, but already at this level of the analysis (i.e. by just focusing on the problems to be 
solved and not on the general problems of method) we pointed out in Subsections 2.4 and 3.4 
that the proper frame is missing and this is not just our own impression. 

There are some calls to theoretical physics from the biology frontier to help systematize 
biology in a rigorous manner (see, e.g., [IE])- We physicists could partially answer these calls 
by, for instance, fitting a given biological problem into a given established physical theory and 
use the logical consistency of the latter. This is an exciting and worth thing to do and when it 
is possible to do that there is plenty of cross-fertilizations between the two fields. In a way our 
proposal of the scar creation-annihilation mechanism for viruses (and the whole "relationship" 
between the Thomson problem, or elastic theory of spherical membranes |41] or the Landau- 
Ginzburg theory of phase transitions [49J on the physics side and virus structure on the biology 
side) is an example. Then the application of physical theories to biological problems might make 
that sector of biology just like that physical theory, hence the frame in that case would be set 
without mayor changes in the methods. This, though, would not solve the problem in general 
but would work only for certain particular features of certain particular biological entities and 
it is not certain at all that this would work all the time (and why should it?). 

We have now identified what is missing: a general theory of physical biology (GTPB), 
namely, the conceptual frame that should always work, for all biological problems, no matter 
the complexity, and that "just" needs to be applied case by case. It must give the correct 
description of the phenomenon and, via mathematical elaborations, must produce predictions. 
This is surely a tremendous task. Let us now see whether this is accomplishable with building 
blocks we have already or whether we need to change something on the physics side as well. 



In 1943 such a change of physics to move the frontier of knowledge from dead to living 
matter was deemed necessary by Schrodinger in his famous lectures delivered at Trinity College 
in Dublin [27] . Two were the main messages there: the carrier of genetic information needs to 
be an aperiodic solid and the thermodynamics of living matter is based on the order from order 
paradigm rather than the order from disorder of dead matter. Both considerations were rooted 
into quantum mechanics. 

The first message was a clear prediction of the structure and role of DNA, discovered only 
ten years later. The second message is still an open question that found some partial answers in 
the work of Prigogine and others who invented Non-Equilibrium Thermodynamics (NET) (see, 
e.g., |5Uj for a textbook chapter). Since living organisms are far from equilibrium open systems, 
NET is surely one of the building blocks we need for GTPB, but it cannot be the end of the 
story because, similar considerations apply to other non-living systems (financial markets, traffic 
jams, weather systems, etc, i.e. what is sometimes called a "complex system") and this is not 
what we are looking for. We need a tailored physics for life, or, in other words, what we might 
call a "principle of distinction" at work exclusively for living matter, and this is still not there. 
Another point is that the role of quantum mechanics as pivotal for life phenomena \27\ [28] needs 
to be fully exploited. 

To build up the GTPB is the most exciting and challenging problem ahead, but it is very 
difficult. For the time being, we can solve the many challenging problems of standard biophysics, 
investigate better the role of quantum physics for life and build up a "GTPV" : a general theory 
of physics for viruses. 
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